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Abstract 

Computational methods for global/local analysis of structures which includes both uncou- 
pled and coupled methods are described. In addition, global/local analysis methodology for 
automatic refinement of incompatible global and local finite element models is developed. Rep- 
resentative structural analysis problems are presented to demonstrate the global/local analysis 
methods. 

Introduction 

The use of composite materials in the design of aircraft structures introduces analytical 
complexity due to the nature of the material systems and their failure modes. Moreover, new 
structural concepts introduce additional modeling complexity due to the details of material and 
geometric discontinuities. Predicting the strength of these structural configurations requires an 
accurate stress analysis capability. However, detailed analysis of complex aircraft structures 
can severely tax even today’s computing environment. Therefore, it is highly desirable to 
utilize detailed modeling only where necessary. Embedding local refinement in a single model of 
the entire structure may lead to highly complex modeling due to the use of transition modeling 
between highly-refined regions and regions with less refinement. In addition, transition modeling 
typically introduces distorted elements into the finite element model which may adversely affect 
the accuracy of the solution. Thus, global/local methods which use different idealizations in 
different structural regions (i.e., global and local regions) have been used to determine the 
detailed response of structures. These global/local methods often reduce modeling complexity, 
however, some of them still require at least some level of transition modeling between the highly- 
refined and less-refined finite element meshes (refs. 1 to 4). For the methods using transition 
modeling, one-to-one nodal correspondence is required on the boundaides of the global and local 
models. Such nodal correspondence will hereafter be referred to as nodal compatibility. 

Global/local methods which do not require such nodal compatibility along common sub- 
domain boundaries provide the desirable feature of modeling flexibility and eliminate the need 
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for transition modeling. These methods fall into one of two categories. In the first category, 
results from a relatively coarse model are interpolated and applied as boundary conditions (i.e.; 
displacements or forces) on an independent detailed local model of the region of interest (see, 
for example, refs. 5 and 6). Thus, methods in this category are considered uncoupled and are 
advantageous when the region of interest is not known prior to or not adequately modeled in 
the initial global analysis. However, typically no provision is made for interaction between the 
global and local models. 

In the second category, the detailed refinement is incorporated within the global model. 
Methods in this category are considered coupled and provide modeling flexibility (i.e., they 
permit independent modeling of the global and local regions) as well as a coupling of the global 
and local analyses (z. e,, they provide the necessary interaction between the global and local 
models). Examples of coupled methods include domain decomposition and coupling techniques 
which focus on the development of such methods for parallel computers (refs. 7 and 8). Other 
methods use some form of multi-point constraints along the global and local boundary (refs. 9 
and 10). 

The purpose of this paper is to describe two methods, one in each of the previously discussed 
categories, that have been developed for global/local analysis. These methods, which include 
both uncoupled (ref. 11) and coupled (ref. 12) global /local methods, have been implemented 
in a general-purpose finite element code, COMET (ref. 13). Both methods eliminate the need 
for transition modeling and nodal compatibility across the common boundary between the lo- 
cal and global models. The uncoupled global/local method uses surface spline interpolation to 
determine the boundary conditions for the incompatible local model. The coupled global/local 
method uses an independent function along the interface between the subdomains to couple the 
incompatible global and local subdomains. In addition, an error detection and control strat- 
egy has been combined with the coupled global/local analysis method to automatically refine 
the independently discretized regions. The description of each of these global/local analysis 
methods including their mathematical formulations is presented first. Then, the methods are 
demonstrated on selected linear, static structural analysis applications. 

Global/Local Analysis Methods 

Global/local analysis procedures which determine the detailed response of structures using 
different idealizations in specific regions are described in this section. Two computational meth- 
ods are described: (1) an uncoupled global/local analysis method and (2) a coupled global/local 
analysis method. In addition, an error detection and control strategy (refs. 14 and 15) is de- 
scribed and is combined with the coupled global/local analysis method in order to automatically 
refine selected independently discretized structural regions. All methods have been implemented 
in a common structural analysis framework (refs. 13 and 17). 

Uncoupled Global/Local Method 

The uncoupled global/local analysis method was developed in detail in reference 11 and 
is briefly described herein. The method employs spline interpolation functions which satisfy 
the linear plate bending equation to determine “boundary conditions” for the local model. 
A schematic which describes the uncoupled global/local analysis method is shown in figure 
1. The global analysis is peiformed first. The interpolation region encompassing the critical 
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region is specified. A surface spline interpolation function is evaluated at every node in the 
interpolation region. The displacement field from the global analysis is used to compute the 
unknown coefficients of the surface spline function. An independent, more refined, local model 
is generated within the previously-defined interpolation region. The interpolated displacement 
field is used to produce a local displacement field which is applied as a “boundary condition” on 
the boundary of the local model. Then, an independent, two-dimensional local finite element 
analysis is performed. 

This method can be used to determine local two-dimensional stress states for specific struc- 
tural regions using independent, refined local models which exploit information provided by 
less-refined global models. The method do.es not require knowledge of the regions requiring 
local detailed stress analysis in advance. However, the analyst may be able to anticipate which 
regions will experience high stresses and the global model can be developed with this knowl- 
edge in mind. In addition, transition modeling is eliminated since independent global and local 
models are used. 

One of the key components of this global/local method is the interpolation of the global 
solution to determine boundary conditions for the local model. The interpolation is performed 
using a surface spline approach. A surface spline is used to interpolate a function of two variables 
and removes the restriction of single variable schemes which require a rectangular array of grid 
points. The derivation of the surface spline interpolation function used herein is based on the 
classical plate bending equation 


D'V‘^io(x,y) = q{x,y). (1) 

where D is the flexural rigidity, q{x,y) is the transverse loading and w{x,y) is the transverse 
deflection of the plate. It ia assumed that q{x,y) = X)"=i l inD ~ ~ " 

is the number of nodes in the interpolation region, pi is the applied point load at the 
node in the region and 6(^x — X{,y — yi') is the two-dimensional Dirac delta function such that 
8(x - Xi,y — yi) — 1, when x = x, and y = yt; and S(x — Xi,y — yi) = 0, otherwise. The 
interpolation function, w(x,y), may be written as 

w(x,y) = ao + aix + 02 y -|- a^x^ -|- a^xy + a^y^ -|- a^x^ -|- ajx^y+ 

a^xi/ -I- -I- ^ F,r] ln(rf ) (2) 

where 

r] = {x - Xif + {y - yif , (3) 

Fi = the coordinates of the node in the interpolation region. The 

logarithmic term, rf In(rf), in equation (2) is the fundamental solution of equation (1) for 
a unit point load at the node in the interpolation region. This type of interpolation was 
first used by Harder and Desmarais (ref. 16) to interpolate pressure loads on idealized wing 
structures. For that application, only three constants oq, <xi, and 02 were used. In reference 11, 
a complete cubic polynomial was used for general structural applications. The n -|- 10 unknown 
spline interpolation coefficients (oq, , 02 ? • • • 1 ^9) Fi) are determined by Evaluating equation (2) 
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at n nodes and solving the resulting equations along with the following additional ten constraint 
equations: 
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These constraint equations are used to prevent equation (2) from becoming unbounded for larp 
values of x and y. The matrix equation is formed by combining equations (2) and ( j, name y 


a p 

0 



where P is an n x 10 submatrix with the i*'* row expressed as 
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Pj = [1 Xi y, 


xi XiXji y, 


Xi 


xi^yi 


Xiyi 


y.^J, 


0 is a 10 X 10 null submatrix, f? is an n x n submatrix with its elements given by 


n 




= { 


-I In('fj). 

0 , 


i ^ j 

i =j 


ij = l,2,...,n 


The subvectors a and F contain the unknown coefficients of the interpolation function, and the 
subvector w contains the generalized displacement field from the interpolation region m t e 
global model. The second set of the equations in equation (5), P F — ® matrix ° 

the constraint equations defined by equations (4). Upon solving equation (5) for the coefficients 
(an fli 02 , . . . , 09 ,F,), equation (2) is evaluated at the I nodes along the global/local interface 
boundary to obtain the interpolated data at the boundary nodes of the local model and is given 

by 

tOg/i{x,,y,) = oo + oix, + 02y. + 03X.2 + aiX,iji + a^yi^ + a^Xi^ + a^x,^yi + 

i = 1 , 2 ,...,/. (6) 


n 


o-sXiyi^ + ogy,^ + ^2, ^ 
i=i 
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If a boundary node of the local model coincides with a node in the interpolation region of the 
global model, this node corresponds to a “knot” in the surface spline and thus, the interpolated 
value is exact. Although equation (6) is derived from the classical plate bending equation, 
it may be used to interpolate each displacem ent component u, u, and w, independently. The 
rotations 0x, and 6z are interpolated (using eq. (6)) from the rotations in the interpolation 
region in the global model instead of differentiating the out-of-plane-deflection, w, to obtain the 
rotations (e.g., io,x). 

Guidelines for the use of the uncoupled global/local method are presented in reference 
11. In this reference, it was shown that the local model boundary should be outside the local 
gradient in order to obtain smooth stresses across the global/local boundary. However, even if 
this requirement is not satisfied, the method provides insight of the local behavior. While the 
method permits simpler modeling than embedding local refinement in a single model analysis, 
it does not take into account any interaction between the global and local models. Thus, an 
alternative method which couples the global and local models together has been developed. 

Coupled Global/Local Analysis Method 

This coupled global/local analysis method (see fig. 2) was developed in detail in reference 12 
and is briefly described herein. In this coupled method, the incompatible global and local models 
are coupled using an independent function along the interface between subdomains. The global 
and local subdomains are discretized by independent finite element meshes. The interface, S , 
(see fig. 3) between the subdomains is discretized independently. The displacement field along 
the interface is expressed in terms of the generalized displacements along the interace. 

The modeling flexibility provided by such a method allows an analyst to incorporate local 
changes in structural design {e.g., hole size and shape) by interchanging different local models 
within the same global model. Moreover, the method may also be advantageous when combined 
with adaptive re^ement strategies. The use of such nodally incompatible models also minimizes 
mesh distortion by avoiding skewed elements otherwise necessary in transition modeling. 

The method described herein may generally be applied to an arbitrary number of indepen- 
dently modeled regions. However, in the following discussion, the mathematical formulation 
will be described in terms of two regions and a multi-segmented interface. A two-dimensional 
domain Q, is modeled as two independently discretized subdomains, Qi and 0 , 2 , as shown in 
figure 3. The interface, S, is modeled as semi-independent straight line segments which conform 
to the configuration but not the discretization of either of the subdomains. Each segment of 
the interface, S, is discretized as a mesh of evenly spaced “pseudo-nodes.” In this paper, the 
following notation conventions are used. The superscript, i, is associated with the mteiface 
degrees of freedom of a subdomain, and the superscript, o, is associated with the other degrees 
of freedom, not on the interface. The subscript, j, denotes the subdomain {e.g., j = 1, 2 for the 
model in fig. 3). Along the A:*'* segment of the interface, S, there are mkj interface nodes for 
subdomain j and n* pseudo- nodes on the interface, 5. Furthermore, it is assumed that each 
node has p degrees of freedom so that there are pm^j degrees of freedom for the interface nodes 
of the subdomain j and pnjt degrees of freedom for the interface, S. 

While three different interface formulations were examined in reference 12, only one is 
presented herein, namely, the hybrid variational formulation. The displacement field, v, along 
segment k of the interface may be written as 
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V - Tq, (7) 

where , T is a p x pnjt matrix containing the interpolating functions, and is a vector of pn^ 
generalized displacements associated with the interface pseudo-nodes. The specific form of T 
depends on the type of function chosen to describe the displacement vector, v, and the number 
of evenly spaced pseudo-nodes, nit, selected along segment k of the interface, S. In reference 12, 
three different functions were examined: linear, quadratic and cubic splines. The cubic spline is 
used in this study. Equation (7) is assumed to be valid along each segment (segments AB and 
BC in fig. 3); at the interface corner (point B in fig. 3), the values from each interface segment 
are constrained to be the same. 

The hybrid variational formulation uses an integral form for the compatibility between the 
interface and the subdomains. The augmented total potential energy is defined as 

n = nn, +nn,+ /a^(v + (8) 

where Aj is a vector of Lagrange multipliers and Uj is the displacement vector along the interface 
for the subdomain j. The total potential energy of subdomain j, Hn . , is given by 
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qjf. 
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; j = i,2 


( 9 ) 


where the generalized displacement vectors q_,, generalized force vectors fj, and the stiffness 
matrices Kj, are partitioned into subvectors and submatrices corresponding to interface nodes 
and other (non-interface) nodes. 

The constraint integrals in equation (8) are added to the functional to enforce, in the 
variational sense, continuity of displacements across the interface. Equation (8) corresponds to 
the “double layer interface” or “frame” method of the hybrid variational principal (ref. 18). In 
this formulation, the Lagrange multipliers. A, are assumed independently at the interface for 
each element, and the displacement vector, v, is assumed independently on the interface, S. 

Assuming that the displacement boundary conditions are satisfied, the stationary condition 
for the total potential energy for arbitrary u in the subdomains, arbitrary v on the interface, 
S, and arbitrary A on the interface parts of the subdomains, results in the following Euler 
equations associated with the interface 


*n = 0 


Aj -l- A 2 — 0 

Uj = v; 


i = l,2 

i = l,2 


on S 


( 10 ) 


These equations are in addition to the usual Euler equations which satisfy the equilibrium 
equations, traction boundary conditions, and the interelement traction continuity. Thus, Xj 
represents the tractions on the interface and the sum of the tractions across the interface is zero 
(i.e., equilibrium is maintained across the interface). The dispacements, Uj, along the interface 
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of the subdomains is expressed in terms of unknown nodal displacements, gj, as Uj = N^q' , and, 
the Lagrange multipliers, Xj, are expressed in terms of unknown coefficients, a^, as Aj = R 
where Nj and Rj are matrices of interpolation functions. The interpolation functions in the 
matrix Rj are taken to be constants for linear elements and linear for quadratic elements. 
Substituting these assumptions into equation (8) and taking the first variation of the potential 
energy with respect to the independent variables, qJjq^jq,, ^md otj, and setting it to zero yields 
the following system of equations: 
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where Mj and Gj axe integrals on the interface defined in terms of Rj, T, and Nj as 

Mj = -jN]Kjds and Gj = J T'^Kjds ; j = l,2 (12) 

Adaptive Global/Local Analysis Procedure 

The adaptive global/local analysis procedure, described by figure 4, is a strategy for auto- 
matically refining selected, independently discretized regions of a structure within the context 
of a linear static analysis. In this method, the enor detection and control strategy presented in 
references 14 and 15 and summarized herein is combined with the coupled global/local method 
previously described. 

Initially, the local and global subdomains are identified and independently discretized. A 
linear stress analysis is performed using the coupled global/local method. Once a displacement 
solution has been obtained for the coupled system, the two models are again treated as inde- 
pendent for post-processing. Refinement indicators, based on stress resultants smoothed only 
within a global or local subdomain, are evaluated, and each subdomain is remeshed as needed. 
If no elements in either model require refining, the solution strategy stops and reports that so- 
lution convergence ha.s been achieved. If either or both models change, a new mesh (or meshes) 
is (are) generated and a new coupled solution is determined. The global and local models are 
parametrically generated using the high-level command language in COMET such that auto- 
matic refinement of the models is accomplished by automatically changing the discretization 
parameters. The strategy continues iteratively until either convergence to a user- specified er- 
ror tolerance is achieved or another riser-specified termination criterion is met (e.g., maximum 
number of mesh iterations or maximum execution time). 

The error detection and control stra.tegy includes the definition of error and refinement 
indicators and the development of a mesh refinement strategy. While several error indicators 
were considered initially, only two proved to be useful in the work presented in references 14 
and 15. The fust is an adaptation of the energy norm of the error; the second is based on the 
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standard deviation in the stress resultants at element nodes. The first error indicator is used in 
the present work and is summarized herein. 

The energy norm of the error as proposed by Zienkiewicz and Zhu (ref. 19) has been 
adapted to plate and shell problems. For the element, the energy norm of the error, He^He, 
is defined by 


IKll|= / (s:-s;rcj‘(s:-siMn.= / asiq'as, 


dO,f 


(13) 


where Sfc is a vector of stress resultant components (where an * denotes smoothed values and 
a ' denotes discrete finite element values); and C is the constitutive matrix that relates stress 
resultants to reference surface strain components. 

As an integrated quantity, the energy norm of the error requires either specific element data 
(i.e., the element shape functions) or an element-independent generalization. This implemen- 
tation provides an element-indepenclent generalization which allows for the decoupling of the 
error analysis from other element-specific functions. The distribution of the difference between 
the smoothed and the discrete values of the stress resultants, AS, over each element is assumed 
to be described by Lagrangian shape functions. This assumption allows for the evaluation of an 
element- independent error indicator and provides a common basis for the evaluation of different 
element types but does not ignore all element -specific features in that discrete element stress 
resultants are calculated using features of the specific element. 

The automatic refinement strategy is based on the local refinement indicators; elements 
needing to be refined are identified using the local refinement indicator (see ref. 19), ik-, defined 
by 

6. = ( 14 ) 

llerlp. 


where 

^111 ^ ntai 

l|u'|( is the energy norm of the current finite element solution over each subdomain, ||e||E is the 
energy norm of the error in each subdomain and is defined as the sum of the elemental values, 
||ejt||p, and Ue/e is the total number of elements. At the start of a specific analysis, the user 
must "specify a value for the percent error in the energy norm which will be acceptable for the 
given analysis; this user-specified tolerance is denoted fjtol- Thus, when ^k is less than unity 
for an element or patch of elements, refinement must occur; when is greater than unity for 
a patch of elements, the patch may be “fused” into a patch containing fewer elements. The 
“fusion” of elements has not been included as part of the error control strategy in this work. 

The error indicator discussed herein requires the definition of a smoothed stress resultant 
field (which is generally availalde in most large-scale structural analysis codes) as an approxi- 
mation to the exact stress resultant field. For simplicity, simple nodal averaging has been used 
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(15) 
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to produce the smoothed solutions for the example cited herein. For all nodes, smoothed stress 
resultants are given by 



JL 



(16) 


where is the smoothed stress component a at node j , rig is the number of elements directly 

connected to node j , and is the discrete stress component a at node j from element L At 
geometric and material discontinuities, the stresses in the same geometric plane or of the same 
material are smoothed separately. It should be noted that any smoothing algorithm may be used 
since obtaining the smoothed solution is a separate part of the error detection strategy. When 
calculating the error and the refinement indicators, it is assumed that a smoothed solution is 
available. 

Automatic mesh refinement may take several forms and need not be fully adaptive. A 
fully adaptive strategy will refine only those elements identified by the refinement indicators as 
elements requiring refinement. Automatic mesh refinement encompasses a range of refinement 
techniques and includes any mesh refinement technique which can be performed with minimal 
user intervention and which can be performed automatically by the structural analysis software. 
The present work (as in references 14 and 15), employs uniform and quasi-uniform refinement 
schemes implemented through a high-level command language facility in COMET. 

Uniform refinement methods may be global but are best employed on parameterized finite 
element models so that only portions of the mesh will be refined and sophisticated data struc- 
tures are not required. This type of refinement is denoted quasi-uniform refinement. Uniform 
and quasi-uniform refinements will tend to over-refine the finite element model. The cost per 
analysis of using such refinement is therefore generally higher than the cost of using a fully 
adaptive refinement technique. Within the context of a general-purpose code, however, a high- 
level refinement technique which is easily implemented, requires no complex data structures, 
and little or no modification to existing software, may well be the technique of choice for many 
analysts. 


Applications 

The effectiveness of the global/local analysis methods described in the previous section 
is demonstrated by determining the response of two structural analysis problems. First, a 
graphite-epoxy blade- stiffened panel demonstrates the uncoupled global/local analysis method 
on a representative aircraft subcomponent with characteristics which often require a global/local 
analysis capability. These characteristics include a discontinuity, eccentric loading, large dis- 
placements, large stress gradients, high in-plane loading and a brittle material system. Sec- 
ond, an orthotropic plate with a central crack demonstrates the coupled global/local analysis 
method. In addition, the plate with a central crack demonstrates the effectiveness of the adap- 
tive global/local analysis procedure. This problem is characterized by a stress singularity at 
the crack tips which presents a challenge for both the global /local method and the adaptive 
procedure. 
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Blade-Stiffened Panel with Discontinuous Stiffener 

An uncoupled global/local analysis of a composite panel with a discontinuous stiffener was 
presented in reference 11 and is summarized herein. A linear stress analysis of a flat, blade- 
stiffened, graphite-epoxy panel loaded in axial compression is performed. The overall panel 
length is 30 in., the overall width is 11.5 in., the stiffener spacing is 4.5 in., the stiffener height 
is 1.4 in., and the cutout radius is 1 in. The cross section of the three blade— shaped stiffeners 
is identical. The loading is uniform axial compression in the x-direction. The loaded ends of 
the panel are clamped and the other sides are free. The material system used for the panel is 
T300/5208 graphite-epoxy unidirectional tapes with a nominal ply thickness of 0.0055 in. The 
panel skin is a 25-ply symmetric laminate ([±45/02/ T 45 /O 3 / ± 45 /O 3 / ± 45/00] and the 
blade stiffeners are 24-ply symmetric laminates ([±45/0io] J. 

Based on the results from a coarse global model analysis, two local regions are considered, 
one near the cutout, and one at the skin-stiffener interface near the edge of the cutout. A four- 
node quadrilateral shell element (ref. 20) was used in the global and local analyses. This element 
is ba^ed on a displacement formulation and includes rotations about the outward normal axis. 
The coarse global model has 2316 active degrees of freedom. The local models of the region 
around the cutout and the skin-stiffener interface have 3456 and 462 active degrees of freedom, 
respectively. The boundary conditions for both local analyses were obtained using data from 
the coarse global model analysis. Two interpolation regions were used for each of the local 
analyses. The first interpolation region, specified in the plane of the panel skin, is used to obtain 
boundary conditions on the local model boundary of the panel skin. The second interpolation 
region, specified in the plane of the stiffener, is used to obtain boundary conditions on the local 
model'boundary of the stiffener. The boundary conditions for the panel skin and the stiffener 
were obtained independently. The displacements and rotations at the skin-stiffener intersection 
on the local model boundaries were taken as those values obtained for the panel skin. 

Contour plots of the longitudinal stress resultant are shown in figure 5 superimposed 
on the deformed global finite element model and undeformed local finite element models. The 
distribution of the longitudinal stress resultant Nj; in the panel skin normalized by the applied 
running load {Nj:)avg (* e., applied load divided by the panel width) is shown in figure 6 as a 
function of the lateral distance from the center of the panel normalized by the radius of the 
cutout. The solution from a very refined global nrodel is used as the reference solution. This 
refined global model lias the same level of refinement in the local regions as the independent 
local models and has 8460 active degrees of freedom. The coarse global model solution provides 
good agreement with the reference solution away from the cutout but underestimates (by 24%) 
the stress concentration at the edge of the cutout. The local model analysis of the region 
around the cutout predicts normalized stress resultants within 1.4% of the refeience solution in 
the neighborhood of the cutout. The local model analysis of the skin— stiffener interface region 
also agrees with the reference solution in predicting a higher gradient than the coarse global 
model analysis. This method may be userl to interrogate any other region of the subcomponent 
as well. 
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Plate with a Central Crack 

An orthotropic plate which has a central crack and which is subjected to uniform tension 
is used to demonstrate the coupled global/local analysis method. This problem is characterized 
by a stress singularity at the crack tips. The remote tensile loading induces mode I crack 
deformations. Accurate solutions are available in the literature (ref. 22) for this configiuation. 

The ratio of the plate half-length to the plate half- width, y, is 2 (see fig. 7). The ratio of 
the crack half-length to the plate half- width, f , ranges from 0.125 to 0.75. It is assumed that 
the material axes (1,2) are along the coordinate axes (x,y). The ratio of the Young’s moduli, 
is 0.1. The ratio of the shear modulus to the longitudinal Young’s modulus , , is 0.8757 

and Poisson’s ratio, i/yx, is 0.21. Using the symmetries in the problem, one quarter of the plate 
is modeled. The local area of interest is first isolated from the larger global model and two 
separate finite element models (shown in fig. 8) are created: one for the global subdomain 
which is very coarse, and one for the local subdomain which is considerably more refined to 
delineate the stress gradient near the crack tip. The finite element nodes of the global and local 
models do not coincide at the global/local interface. The local model dimension is defined by 
the parameter, H, shown in figure 7. A nine-node Assumed Natmal-Coordinate Strain element 
21 was used in this problem. This element has five degrees of freedom at each node (i.e., three 
displacements, two bending rotations and no drilling degrees of freedom) and uses a strain 
approximation to calculate the element stiffness matrix (which is equivalent to a a stiffness 
matrix calculated using a selectively directionally reduced order of integration). A linear static 
analysis is performed using the coupled global/local method to join these two incompatible finite 
element models. 

The stress intensity factors, A", were calculated from the relationship 


K = 



(17) 


where G is the strain-energy release rate and c is a constant defined in terms of the material 
' properties. The strain-energy release rate, G, is calculated by the virtual crack-closure technique 
and for quadratic elements (ref. 23) is defined as 


G = -■^[FjVj^2 + Fj-iVj^i\ (18) 

where A is the element size at the crack tip, Fj and Fj-i are internal forces at the crack tip and 
one node point ahead of the crack tip, respectively, and Vj^\ and Uj+2 are the displacements at 
the first and second node points behind the crack tip (*.e., on the crack face, see fig. 9). Note 
that no special element (e.y., quarter-point or other singularity element) is used at the crack 
tip. The value of c for orthotropic materials is given in reference 24 as 

^ 11^22 \ ^ 

2 ) 
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A 22 2Ai 2 + Aee 
+ 2A,i 


( 19 ) 


where 
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The stress intensity factors, K , normalized by dQy/ifa are shown in figure 10 as a function 
of the crack half-length to plate half- width ratio, f. The solid curve in the figure represents 
the solution by Bowie (ref. 22). The open circles represent the stress intensity factors obtained 
from the finite element analyses using the coupled global/local method. The calculated stress 
intensity factors are in excellent agreement with the solution by Bowie (ref. 22). 

An adaptive analysis is performed on the plate with a central crack by again isolating the 
local area of interest from the larger global model and creating two separate finite element 
models: one for the global subdomain which is initially very coarse, and one for the local 
subdomain which is also initially coarse. A linear static analysis is again performed using the 
coupled global/local method to join the two incompatible finite element models. Refinement 
indicators based on the energy norm of error, with a 5% error tolerance, are evaluated on Model 
1 (shown on the upper left of fig. 11 with an enlarged local subdomain). These indicators 
suggest that refinement of both the global and local models is required (darker shaded elements 
on the figure indicate that refinement is needed). New models are generated automatically 
for both global and locaT subdomains (denoted Model 2 in the center of the upper portion of 
the figure), an analysis is performed, and refinement indicators are evaluated. In subsequent 
analyses, refinement is required only in the local subdomain of Models 2 and 3 as shown by the 
refinement indicators in the figure. A new local model is automatically generated, the global 
model is allowed to remain the same, and a fourth and final linear analysis is performed on 
Model 4 (lower left of the figure). Although refinement indicators again suggest that refinement 
is needed, the solution procedure is terminated. As the mesh becomes more refined in the 
neighborhood of the crack tip, the solution is driven to an increasingly better representation of 
the singularity, thus the refinement will continue for this problem until user- terminated. The 
transverse stress resultant, Ny, along the width of the panel at the centerline (along the x 
axis shown in fig. 7) is shown in figure 12 for each of the four models. The sequence of four 
models generates increasing values of peak stress as the solution is driven to a better numerical 
representation of infinity at the crack tip. The stress intensity factors predicted using the four 
models are shown in figure 13 as a function of the number of degrees of freedom in each model. 
The stress intensity factor increases and approaches the reference solution by Bowie 22 as the 
number of degrees of freedom is increased. The stress intensity factor predicted by Model 4 
(1605 degrees of freedom) is within 0.1% of the Bowie’s solution. 

Concluding Remarks 

In summary, two computational methods for global/local analysis of structural components 
have been described. The methods are uncoupled and coupled global/local analysis methods. 
Also, methodology for automatic refinement of incompatible global and local models has been 
developed. The effectiveness of the methods has been demonstrated using selected linear static 
structural analysis applications. 

The uncoupled global/local analysis method presented herein reduces the analyst’s compu- 
tational effort and removes the restriction of having to know the regions requiring local detailed 
stress analysis in advance. The coupled global/local analysis method provides interaction be- 
tween the global and local models. Both global/local methods described herein eliminate the 
need for potentially complex transition modeling from highly-refined regions to regions with less 
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refinement and thus reduce the modeling effort. The global/local analysis capabilities provide 
efficient strategies for predicting local detailed stress states using relatively coarse finite element 
models in the far-field. The error detection and control strategy combined with the coupled 
global/local analysis method has been shown to be an effective procedure for improving the 
accuracy of global and local results. The combined procedure provides a modeling technique 
for automatically refining incompatible global and local models. 

The global/local analysis methods described and demonstrated herein, provide techniques 
for predicting local detailed structural response. In addition, the use of these methods should 
reduce the analyst’s modeling effort and thus, should favorably impact the design process for 
structural components. 
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Figure 1. Description of Uncoupled Global/Local Analysis Method. 



Figure 2. Description of Coupled Global/Local Analysis Method. 
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Figure 4. Description of Adaptive Global/Local Analysis Method. 
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Figure 5. Contour Plots of Longitudinal In-plane Stress Resultant, N^. 
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Figur6 6. Longitudinal In-pla.ne Stress Resultant, Nxy Distributions. 
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Figure 7. Orthotropic Plate with a Central Crack. 


y ♦ 


H H 


m 





H 


Figure 8. Finite Element Models for One Quadr ant of Plate with Central Crack. 
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